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Here we analyze the behavior of dynamically coupled maps, based on those introduced in a series 
of papers by Ito & Kaneko (Phys. Rev. Lett., 88, 2002 & Phys. Rev. E, 67, 2003). We show how the 
microscopic coupling mechanism changes the behavior of the system both by affecting the stability 
of fixed points and through a more subtle effect in the crossover behavior between different regions 
of the parameter space. This makes it necessary to choose very carefully the exact manner in which 
one couples maps if they are to be used as a general model of composite systems. 



I. Introduction 

It is generally accepted that the specific details of mi- 
croscopic couplings are irrelevant for the collective be- 
havior of globally coupled maps (GCMs) (see [3] and ref- 
erences therein). There are a multitude of arbitrarily 
different ways in which maps may be coupled together, 
with none being clearly and generically better. We show 
here that the exact mechanism of coupling can lead to 
distinct behaviors of the system. Thus making it neces- 
sary to choose very carefully the exact manner in which 
one couples maps if they are to be used as a general model 
for composite systems. 

Synchronization in systems of coupled maps has been 
widely observed [TJ [21 |4j EJ [7] . They are therefore of 
immense interest as models for understanding synchro- 
nization phenomena as observed in the real-world such 
as fireflies [5], heart pacemaker cells [3] etc. However, 
many of the models to date have their connections pre- 
scribed a priori and as such they are fixed in time. This 
is in contrast to real-world systems where both the units 
and the connections are dynamic elements. It is therefore 
desirable to extend these coupled systems so as to allow 
the interactions to evolve along with the node dynamics. 

Models of interacting iterative maps fall into two gen- 
eral classes: The first consists of those systems where the 
maps are coupled prior to the nonlinear transformation, 
as in 

< +1 = /((l-e)4 + eE /l (4)) C 1 ) 
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which shall be referred to as internal coupling (see for 
example [HE]). 

Alternatively, the nonlinear transformation may be ap- 
plied prior to the coupling, as in 
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which shall be referred to as external coupling (see for 
example d El [TD]). 

In the case when h(x) — f(x) for external coupling and 
h(x) — x for internal coupling, the systems show equiva- 
lent dynamics. The behavior of most interest is that for 
sufficiently high coupling strength, the node-states syn- 
chronize; often providing an example of the surprising 
phenomenon that is synchronized chaos. We show here 
that it is possible to switch between this synchronized 
chaotic behavior to a synchronized stationary state de- 
pending on the way the coupling is implemented, namely 
the functional form of h(x). 

II. Model Details 

The models studied here are based on those introduced 
by Ito & Kaneko in [TJ |2]. When external coupling is 
employed, the system is described by 
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The functional form of g was chosen so as to employ 
Hebbian dynamics |llj . although the results shown here 
are qualitatively unchanged if anti-Hebbian dynamics are 
employed by using g(xt,x 3 n ) = \x % n — 41 instead of Eq. 
([5]) (compare Figures [I] and [6]). 6 is a constant that 
governs the plasticity of the network and is set to 0.1 
throughout the present work. 

The equivalent model with internal coupling is defined 
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by the same equations, only with Eq. ([3| replaced by 
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As with many GCM studies, we use the logistic map, 
f(x) = ax(l — x) as the underlying nonlinear dynamics. 
This is one of the simplest functions that can display such 
distinct behavior, from stationary-state, through peri- 
odic to chaotic. There has also been recent evidence to 
suggest that functions with the form of the logistic map 
may indeed be of direct relevance in neuroscience |12j . 
We will consider two different forms of h{x), 

h{x) = x (7) 

h(x) = f(x) = ax{\ ~ x) (8) 

When we use h{x) = f(x) in the externally coupled sys- 
tem, the model is similar to many externally coupled sys- 
tems in the literature. Yet, this corresponds to interac- 
tions between the units occurring instantaneously which 
is obviously not possible in real systems. We therefore 
introduce a degree of causality to the model by simply 
using h(x) = x instead. 

Each system has two basic parameters governing the 
dynamics: the nonlinear parameter a and the coupling 
strength c. This a-c parameter space is known to have 
a series of regions of qualitatively different behavior (see 
[TJ [2] for a detailed discussion of these) . These regions 
are present in both the internally and externally coupled 
systems studied, although the boundaries occur at dif- 
ferent actual values and with different characteristics as 
will be shown in Section IIII CI 

III. Results 

As was found for systems with fixed couplings [TUHT^] . 
we find distinct differences in the dynamics of the three 
systems: The system with external coupling and h(x) = 
x is found to enter a stationary-state after an initial tran- 
sient time, see Fig. [I] This model could be thought of as 
representing systems where interactions occur at finite 
speed and there is therefore a time lag involved in the 
updating. The value of the node-states in the stationary 
state corresponds to the fixed point of the logistic map 

X = 2=1. 
a 

In contrast, both the system with internal coupling 
and that with external coupling but h(x) = f(x) are 
never found to evolve into this state. The behavior in 
the coherent region appears to vary at random, see Fig. 

m 

This distinction is not immediately obvious since the 
fixed point of the logistic map, x — ^=1 is a fixed point of 
both the externally and internally coupled map systems. 
It must therefore be the stability of this fixed point that 
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FIG. 1: Evolution of the node-states (x 1 ) for the external 
coupling, h(x) = x and a — 3.97, c = 0.5, N = 100. The 
time-evolution of each node-state, x % n , is represented by a line 
on the graph. Thus giving N lines on each graph. However, 
when the nodes are synchronized these all lie on top of one 
another so only one line may be seen. The nodes clearly 
synchronize after ~ 1000 timesteps. 



is crucially dependant on the coupling mechanism em- 
ployed. 

In addition, the boundaries between the different 
phases of parameter space occur at slightly different val- 
ues of the parameters and with different characteristics. 
This can be seen by comparing Figs. [7] [8] and [9] Once 
again there is a stark difference between the behavior of 
the externally coupled system with h(x) — x and the 
others. The trend with increasing system size for the ex- 
ternal system and h(x) = x is for the boundary between 
ordered and coherent behavior to occur at lower coupling 
strength for larger systems. This is the opposite of what 
happens with increasing system size for both the inter- 
nally coupled system and the externally coupled system 
when h{x) — f(x) where the boundary between ordered 
and coherent behavior shifts to higher coupling strength 
for larger systems. 

A. Fixed Point Stability: Analytical Results 

The stability of a fixed point of a dynamical system 
can be analyzed through a simple perturbation analysis. 
Even to first order we can see where the distinction be- 
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constructively. This cancelation of one another is what 
leads to the reduction in 8 l n and therefore the stability of 
the fixed point. 

We have analyzed these terms through numerical simu- 
lations and found the results (not shown here) to confirm 
this conjecture of the terms canceling one another for the 
external system with h(x) = x. For this system, the two 
terms were found to be of comparable amplitude and op- 
posite sign and thus the fixed point of the system is stable 
to perturbations. 



B. Fixed Point Stability: Numerical Results 



FIG. 2: Same as Fig. [T] but for the model with internal 
coupling and h(x) — f(x) for a = 3.97, c = 0.6, N = 100. 



tween the externally coupled systems with h(x) = x and 
h(x) — f(x) comes from. For a fixed point to be stable, 
any perturbation, S l n = (x l n — x) must decay with time. 

For the internally coupled system with h{x) = x, to 
first order we have 
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For the externally coupled system with h(x) = f(x), to 
first order we have 



J n+1 



h(x)=f{x) 



(2-a)(l-c)4 + (2-o)c(X>«# 

(10) 



For the externally coupled system with h(x) — x, to first 
order we have 



h(x)=x 



(2-a)(l-cK + c(^<^ 



(11) 



Consider Eq. ( |T0| ): The pre- factor (2 — a) to the second 
term plays a vital role in the region we are interested in, 
a > 2. Here, (2 — a) < so the first and second terms 
will be of the same sign and therefore add constructively 
causing 8 l n to oscillate about the fixed point. Namely, 
that if S l n < then S l n+1 > 0. Whereas if 5 % n > then 

S l n+i < 



In contrast, for the system with h(x) = x, Eq. (Ill 
there is no negative pre-factor to the second term. Thus 
the two terms will be of opposite sign and so add de- 



FIG. 3: A bifurcation diagram for the system with external 
coupling and h(x)—x. The plot shows the behavior of all 
nodes for 100 time steps at various c-values for fixed a = 3.97. 
For each time step, each x l n is plotted. For higher coupling 
values when the system is synchronized, these will all coincide 
and there will therefore appear to be less points plotted. 

Through numerical simulations of the three systems we 
have been able to characterize the behavior displayed by 
each system. Through plotting bifurcation diagrams, we 
can see how this behavior is dependant upon the coupling 
strength parameter, c. 

Fig. [3] shows us that the externally coupled system 
with h(x) = x is in fact stable on the fixed point for a 
range of coupling strengths (c-values). This is in contrast 
to the other systems whose bifurcation diagrams show 
that the systems do not evolve onto a fixed point at any 
coupling strength, see Figs. [4] and [5] 

When the system is reduced to the standard GCM with 
all-to-all couplings that are constant in time, these dis- 
tinguishing behaviors of the systems are effectively un- 
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FIG. 4: A bifurcation diagram for the system with external 
coupling and h(x)=f(x). The plot shows the behavior of all 
nodes for 100 time steps at various c-values for fixed o = 3.97. 
The points become markedly less dense after c ~ 0.5 because 
the system is synchronized in this region. Thus, each time- 
step results in only one point on the graph as opposed to 
100 (one for each node) as is the case when the system is 
unsynchronized . 




FIG. 6: The same as Figure [4] but employing Anti-Hebbian 
dynamics in the evolution of the connection strengths. The 
plot shows the behavior of all nodes for 70 time steps. 



changed. Namely, that the externally coupled system 
with h(x) = x has a stable fixed point whereas the other 
systems do not. This result should be expected since the 
stability as explained through linear perturbation analy- 
sis is not a result of the network topology. 



C. Phase Boundaries 



1.0 

___ o.s 

CD 0.6 
03 

? 0.4 


O 0.2 







"I. 
[' 

"'A, 


! i; =.!>': 

: ?" ;,■ 
I • - a ' 


S , « I , ; ■ 










1 fjri 








.. 




';!->' 




f.-i-' " 


P ! 
















ZhlxZ Sl„l|iM r",t"/ 



0.2 0.4 0.6 0.8 1.0 

C 



FIG. 5: A bifurcation diagram for the system with internal 
coupling. The plot shows the behavior of all nodes for 100 
time steps at various c-values for fixed a — 3.97. 



As was shown by Ito & Kaneko in [5] the a-c parameter 
space consists of distinct regions which are qualitatively 
different. The regions are essentially those first described 
in |4J ; namely that for sufficiently high coupling strength 
there is a region of a-c parameter space where all nodes 
synchronize, displaying coherent behavior. For lower cou- 
pling strength, this global synchronization is lost and the 
system splits into distinct clusters, this region of a-c pa- 
rameter space is known as the ordered region. For yet 
lower values of the coupling strength there is a disordered 
region, where there is no synchronous behavior between 
any pair of nodes. 

We have found that as one crosses the boundary from 
ordered to coherent regions, the different systems show 
qualitatively different transitions. This can be seen by 
looking at the behavior of a 2 as we cross the boundary 
by changing c at fixed a-value. Where, a 2 is as defined 
in [14]. Namely, 

- 2 = iy(EK-<^>] 2 ) t (12) 
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where < x n > denotes the average over all node-states 
x 1 at timestep n and (•••), denotes the average over all 
timesteps. 




FIG. 7: Here we see the transition to synchronization and 
how it changes with the system size for the internally coupled 
system. For ease of comparison, we re-scale the data and plot 
a 2 versus the coupling strength c. This data is the average 
over 1000 random initial conditions. Here we show the tran- 
sition for different system sizes, N = 10 (□), N = 100 (O), 
N = 500 (x), N = 1000 (+) and N = 2000 (A). 



Figs. [7] [8] & [9] show that there is a fundamental differ- 
ence in the crossover from ordered (distinct synchronized 
clusters) to coherent (one giant synchronized cluster) re- 
gions between the three systems. The externally cou- 
pled system with h{x) = f(x) and the internally coupled 
system have qualitatively similar behavior in er 2 ; with 
increasing system size the transition moves to higher c- 
values and the sharpness increases. There is however one 
local defining feature for the internally coupled system, 
at c = 0.5: there is a localized decrease in a 2 . In con- 
trast to either of these two systems, for the externally 
coupled system with h(x) — x, the transition moves to 
lower c-values for larger system sizes. If these results 
are extrapolated to the limit of infinite system size, this 
suggests that only the externally coupled system with 
h(x) = x would be able to synchronize globally. 

One feature that is common to all three systems how- 
ever is the fact that the nature of the transition from one 
phase to another (ordered to coherent) is fundamentally 
different to other types of phase transition. A hallmark 
of phase transitions in equilibrium statistical mechanics 
is a discontinuity in the order parameter, in the limit of 
an infinite system size. Numerically, this can be seen 
by a sharpening of the transition with increasing system 
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FIG. 9: Same as Figure [7] but for the externally coupled sys- 
tem with h(x) — x. 
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size. However, here the sharpening as seen in Figs. [7]& 
[8] possesses a subtle distinction since it is a probabilis- 
tic measure of the size of the basin of attraction for the 
coherent state. 

For a given point in the a-c parameter space and set of 
initial conditions, the system will evolve into the coherent 
state say. For an arbitrarily small change in the initial 
conditions we see that the system is attracted to the or- 
dered state of several clusters. This apparent riddling of 
the basin of attraction of the coherent region is ordinarily 
not observed in standard phase transitions and gives rise 
to a broad and structured boundary between the coherent 
and ordered phases in the a-c parameter space. Indeed, 
such a riddling has been found in other coupled map sys- 
tems [IS]. It has also been long since conjectured that 
riddled basins of attraction could be extremely common 
in natural systems |16j . In such a scenario, the param- 
eter space phase diagrams can be misleading since they 
give no indication of this broad parameter region that is 
neither wholly in one phase nor the other. 

IV. Conclusions 

We have shown that the microscopic details of dynamic 
couplings have a dramatic effect on the macroscopic dy- 
namics displayed by systems of globally coupled maps. 
Systems with asynchronous coupling are, in some cases, 
able to stabilize a fixed point of the underlying map. This 
stabilizing effect of coupling shown here is the same as 
was previously shown for coupled map lattices with fixed 



un-weighted coupling (TUl H3]- Since the stabilizing ef- 
fect is found when there is a finite speed at which the 
interaction between nodes is transmitted, this does call 
into question the applicability of other coupled systems 
as prototype models for synchronized chaos of real world 
systems. 

We must also question what, if any, significance may be 
attached to the form of coupling that must be employed 
in order to gain synchronized chaos, namely the internal 
coupling or synchronously updated external coupling. If 
these models are to be viewed as prototype models of 
real world synchronization phenomena, we may use this 
result to gain insight into the mechanisms governing real 
systems. 

Further, we have shown that this stabilizing of fixed 
points is not the only macroscopic difference between the 
systems. The crossover from one region of parameter 
space to another has been shown to be highly variable 
between the different systems studied. Although all sys- 
tems show non-trivial mechanisms as one crosses from 
one phase to another. The boundary maintains a finite 
width with a great deal of structure to it in all three 
systems studied. 
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